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We point out that events associated with large Et direct photons in high-energy heavy-ion collisions 
can be used to study jet energy loss in dense matter. In such events, the pr spectrum of charged 
hadrons from jet fragmentation in the opposite direction of the tagged photon is estimated to be well 
above the background which can be reliably subtracted at moderately large px- We demonstrate 
that comparison between the extracted fragmentation function in AA and pp collisions can be used 
to determine the jet energy loss and the interaction mean-free-path in the dense matter produced 
in high-energy heavy-ion collisions. 
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Large transverse momentum jets, among many other 
hard processes in high-energy heavy-ion colhsions, have 
been proposed as effective probes of the transient dense 
matter. For example, an enhanced acoplanarity and en- 
ergy imbalance of two back-to-back jets can be used to 
study multiple scatterings of a parton inside a dense 
medium |l[]. Study of large px jets can also probe their 
energy loss due to inelastic scatterings inside a dense 
matter or a quark-gluon plasma Q . Because of the enor- 
mous background in high-energy heavy-ion collisions, the 
conventional calorimetric study cannot measure the jet 
energy to such an accuracy as required to determine 
the energy loss. Alternately, single-particle inclusive pr- 
spectrum has been shown to be sensitive to the jet en- 
ergy loss [^. Since the single-particle spectrum is a 
convolution of the jet production cross section and the 
jet fragmentation functions, the suppression of produced 
hadrons at a fixed pT results from jet quenching with a 
wide range of initial transverse energies, thus making it 
difficult to measure directly the modification of jet frag- 
mentation for a given transverse energy. 

In this Letter, we propose to study the jet quench- 
ing in high-energy heavy-ion collisions by measuring 
the pt distribution of charged hadrons in the oppo- 
site direction of a tagged direct photon. A direct 
photon is produced by quark-antiquark annihilation or 
quark(antiquark)-gluon Compton scatterings in which a 
gluon or quark(antiquark) jet is also produced in the op- 
posite direction of the photon. By tagging a direct pho- 
ton with a given transverse energy E^, one can avoid the 
uncertainties associated with the jet production cross sec- 
tion. One can also determine the initial transverse energy 
of the produced jet, ~ E'^, from momentum conser- 
vation, modulo calculable corrections from initial state 
radiations. At collider energies and sufficiently large E^, 



the Cronin effect due to multiple scatterings during the 
initial interaction stage is also negligible jj]. We shall 
use perturbative QCD to show that the pr spectrum of 
charged hadrons with moderate px in the backward di- 
rection of a direct photon is a very good approximation of 
the jet fragmentation function which can thus be reliably 
extracted. We shall also study the sensitivity of the mod- 
ification of the jet fragmentation functions in heavy-ion 
collisions to the energy loss of jets and the jet interaction 
mean-free-path inside a dense matter. 

Let us consider events with a direct photon in the cen- 
tral rapidity region, \y\ < Ay/2, Ay = 1. For suffi- 
ciently large i?J of the photon, the rapidity distribution 
of the associated jet is also centered around zero rapid- 
ity with a comparable width. If the azimuthal angle of 
the photon is 0-,. and (p^ ~ (p^ -\- n, most of the hadrons 
from the jet fragmentation will fall into the kinematic 
region, {\y\ < Ay/2,\(f>— < A(j)/2), where one can 
take A(p — 2 according to the jet profile as measured in 
high-energy collisions Given the jet fragmentation 
functions Dh/a{z), with z the fractional momenta of the 
hadrons, one can calculate the differential distribution 
of hadrons from the jet fragmentation in the kinematical 
region (Ay, A(/i), 



dycPpT 



Dn/a{PTlET)C{Ay,A4>) 



a.h 



PtEj 



AyA4> 



(1) 



where C{Ay, A<j>) = J\y\<^y/2 dy 
(j)~f) is an overall factor and f{y, 4> 



is the hadron profile 
aroimd the jet axis. The summation is over both jet (a) 
and hadron species (/i), and ra{E^) is the fractional pro- 
duction cross section of a-type jet associated with the 
direct photon. We define D^{z) — '}Zah'''o.iE'^)Eh/a{z) 
as the inclusive fragmentation hmction. C{Ay,A(j)) is 
the acceptance factor for finding the jet fragments in 
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the given kinematic range. We find C(Ay, A(/)) « 0.5 
at ^/s = 200 GeV, independent of the photon energy 
using HIJING D Monte Carlo simulations for the given 
kinematic cuts. For a fixed E^, the jet Ex has a smearing 
around E^ caused by initial state radiations. One should 
therefore average Eq. (|^) over such a smearing. The re- 
sultant spectrum is very well approximated by Eq. (|l|) 
with Et — E^ [Ql, as will be shown by comparison with 
explicit HIJING Monte Carlo simulations. 

To calculate the background for the photon-tagged jet 
fragmentation from particle production in a normal cen- 
tral nucleus- nucleus collision, one convolutes the frag- 
mentation functions with the jet cross sections |q], 



before they escape and fragment into hadrons. There 
have been many studies on the energy loss of a propa- 
gating parton inside a medium. It is believed that ra- 
diative energy loss dominates even when the Landau- 
Pomeranchuk-Migdal suppression is taken into account 
||l3|,^. While a dynamical study of the jet propaga- 
tion and the modification of the hadronization is more 
desirable, we will use a phenomenological model here to 
demonstrate how sensitive our proposed measurement in 
the photon-tagged events to the interactions and the av- 
erage energy loss suffered by a parton in a dense medium. 

We restrict ourselves to the central rapidity region so 
that a parton will only propagate in the transverse 
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= Xriey/Xa + e y/xb)/2, Xf„nin = 

XaXre'"/ (2xa - XTe^), Xamin = xte" / (2 - XTe~y), and 
XT = 2pt/\/s. The K 2 factor accounts for higher 
order corrections The parton distribution density 
in a nucleus, fa/A{x,r) = tA{r)Sa/A{x,r)fa/N{x), is as- 
sumed to be factorizable into the nuclear thickness func- 
tion tA^j) (with normalization J d'^rtAir) — A), parton 
distribution in a nucleon fa/Nix) and the parton shad- 
owing factor Sa/Aix,r) which we take the parametriza- 
tion used in HIJING model §. In our notation, the 
scale dependences of the parton distributions fa/N{x, Q^) 
and the fragmentation functions Di^/i^{z, Q^) are implicit, 
which we take to be Q = E^. 

Jet fragmentation functions have been studied exten- 
sively in pp, ep and e+e~ experiments [|lO|. We will use Di^^g^{z, AL, Q^) 
the parametrizations of both z and dependence of the 
most recent analysis []lT| for the unmodified fragmenta- 
tion functions D'^^^(z), in which only pions and kaons 
are included. We will use the MRS D— ' parametrization 
of the parton distributions [O. The resultant single- 
particle pt spectra from Eq. (0) for pp and pp collisions at 
different energies agree well with the experimental data 
at moderate pr > 2 GeV/c @] where particle produc- 
tion from soft processes is expected to be small. Shown 
in Fig. 1 are the differential px distributions from the 
fragmentation of a photon-tagged jet with Ej, =15, 20 
GeV and the underlying background of normal central 
Au + Au collisions at ^/s — 200 GeV. The points are 
HIJING simulations of lOK events and solid lines are 
numerical results of Eqs. (||) and (||), in both cases no 
medium effects have been considered in the fragmenta- 
tion functions. The background in pp collisions is about 
1200 times smaller than Au + Au. 

In heavy-ion collisions, produced partons will experi- 
ence secondary scatterings and induced radiation which 
will drive the system toward equilibrium. As a result, 
large momentum partons will lose part of their energy 



direction in a cylindrical system. The parton will not 
hadronize inside a deconfined quark-gluon plasma. In 
a hadronic medium, we assume that the fragmentation 
functions can be approximated by their forms in vacuum. 
We only study the effects of radiative energy loss. Given 
the inelastic scattering mean- free-path, , the probabil- 
ity for a parton to scatter n times within a distance AL 
before it escapes the system is assumed to be 
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If we assume the average energy loss per scattering suf- 
fered by the parton is Ca, the modified fragmentation 
functions can be approximated as, 
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where z° = z/{l — nea/Er), z'^ — zEr/ca and C% — 
J2n=o -Pti{n) . We hmit the number of inelastic scatter- 
ings to iV = Et / ta by energy conservation. For large 
values of N , the average number of scatterings within a 
distance AL is approximately {ua) ~ AL/Aq. The first 
term corresponds to the fragmentation of the leading par- 
tons with reduced energy Et — nea and the second term 
comes from the emitted gluons each having energy on 
the average. For simplification, we have neglected the 
fiuctuation in the energy carried by each emitted gluon 
and its possible rescatterings. The inelastic scatterings 
suffered by the leading parton are normally not hard. 
Therefore, we also assume the scales in the fragmentation 
functions of the emitted gluons are given by the initial 
value Qq — 2.0 GeV'^. Since the emitted gluons will only 
contribute to hadrons with very small fractional energy, 
the final modified fragmentation function in the moder- 
ately large z region is not sensitive to the actual radiation 
spectrum and the scale dependence of the fragmentation. 
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FIG. 1. The differential pr spectrum of charged particles 
from the fragmentation of a photon-tagged jet with =15, 
20 GeV and the underlying background in central Au + Au 
collisions at ^/s = 200 GeV. The direct photon is restricted to 
\y\ < Ay/2 — 0.5. Charged particles are limited to the same 
rapidity range and in the opposite direction of the photon, 
10 — 0^ — 7r| < A0/2 = 1.0. Solid lines are perturbative calcu- 
lations and points are HIJING simulations of lOK events. The 
dashed lines are calculations with jet energy loss, dEq/dx — 1 
GeV/fm and the mean- free-path A, = 1 fm. 

Since the jet production rate is proportional to the 
number of binary nucleon-nucleon collisions, the averaged 
inclusive fragmentation function of a photon-tagged jet 
in a central nucleus-nucleus collision is 



DlA^) = J ^J^^r,(£;J)i?,/„(z,AL) , (5) 

ah 



where Taa{0) = J (Prt\{r) is the overlap function of 
AA collisions at zero impact-parameter. Neglecting the 
transverse expansion, AL(r, (f> — (p^) only depends on the 
jet production position {r,(j)). Using Eq. (|4) in Eq. (||), 
we can calculate the single-particle inclusive pt spectrum 
of normal central AA collisions taking into account jet 
quenching. 

In principle, Sa and Aa are related to each other in a dy- 
namical model [p"3|jr^ ] . Phenonienologically, we can treat 
them as independent parameters. Alternatively, we will 
vary Aq and dEa/dx = ea/^a in our calculations. The 
dashed lines in Fig. 1 are calculated with the modified 
fragmentation functions, with dEq/dx = 1 GeV/fm and 
Xq — 1 fm. Wc have assumed that the mean-free-path 
of a gluon is half and the energy loss is twice that of a 
quark. During the parton propagation, multiple scatter- 
ings can also change the direction of the parton result- 



ing in a sizable acoplanarity. Such an acoplanarity due 
to multiple scatterings is probably small as compared to 
that caused by initial state radiations for a large E^ pho- 
ton. Thus, we assume the acceptance factor C{Ay, A0) 
to be the same as in pp collisions. One observes that 
there is significant suppression of large pr particles both 
from the background and jet fragmentation in the oppo- 
site direction of a tagged photon due to jet quenching. 
Since the number of particles at large pt > 4 GeV/c 
from the underlying background is substantially smaller 
than from the tagged jet fragmentation with and without 
jet quenching, one can accurately measure the jet frag- 
mentation function from the px distribution of charged 
particles in the opposite direction of the tagged photons, 
given enough number of events. Once the background is 
subtracted, one can push the limit to even smaller pt > 2 
GeV/c, which corresponds to z ^ 0.1. One can then com- 
pare the fragmentation functions measured in pp, pA or 
peripheral AA with central AA collisions to obtain the 
modification due to jet quenching. 
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FIG. 2. Ratio of the inclusive fragmentation function of 
a photon-tagged jet with and without energy loss in central 
Au + Au collisions for a fixed dEq/dx — 1 GeV/fm. 

To study the sensitivity of the modified inclusive frag- 
mentation function to the energy loss, and the interac- 
tion mean-free-path. A, we plot in Fig. 2 the ratio of the 
fragmentation functions with and without energy loss for 
central Au + Au collisions. There is enhancement of soft 
particle production due to induced emissions, but only 
at very small values of z. The fragmentation function is 
suppressed for large range of z due to energy loss. For 
fixed dEq/dx — 1 GeV/fm, the suppression is delayed to 
larger values of z for larger jet energies. The most opti- 
mal situation is when the average total energy loss (AEx) 
is comparable to the initial jet energy so that substan- 
tial suppression happens at moderate values of z. From 
Figs. 2 and 1, we can see that there is such a window of 
opportunity between = 10 and 20 GeV at y/s = 200 
GeV where the background is small. 

For large values of z > 0.9, particles from the leading 
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jets, which have suffered at least one inelastic scattering, 
are completely suppressed. The remaining contribution 
comes from only those jets that escape the system with- 
out a single scattering. The suppression factor is given by 
(exp(— AL/Aa)), independent of jet energy Et and the 
energy loss dEa/dx. Therefore, one can determine the jet 
interaction mean-free-path by measuring the suppression 
factor of the jet fragmentation function at large z > 0.9. 
For intermediate values of z ^ 0.2-0.5, particles from the 
leading partons with reduced energy dominate as far as 
Ea <IC Et, a situation we will refer to as the "soft emis- 
sion" scenario. Since the average total energy loss by the 
leading parton is (AEto) = {na)ia = {AL)dEa/dx, the 
suppression factor should scale with dEa/dx, depending 
very weakly on the mean-free-path. Shown in Fig. 3 are 
the suppression factor at z = 0.3 as a function of dEq/dx 
for three different values of the mean-free-path. We see 
that for the soft emission scenario, the suppression fac- 
tor scales and decreases almost linearly with dEq/dx. At 
large values of dEq/dx and A^, the average total energy 
loss becomes comparable or equal to the initial jet energy 
Et- In this "hard emission" scenario, particle production 
from the emitted gluons and contributions from those jet 
partons which escape the system without any induced ra- 
diation become important. This is why the suppression 
factor saturates at larger values of dEq/dx, especially for 
large values of A^. Since the mean-free- path can be de- 
termined from the measured suppression factor at large 
z > 0.9 which is independent of dEq/dx, additional mea- 
surements of the suppression at intermediate z = 0.2-0.4 
will enable one to extract the energy loss. 
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FIG. 3. Ratio of the inclusive fragmentation function of 
a photon-tagged jet with and without energy loss in central 
Au -\- Au collisions at 2 = 0.3 as a function of dEq/dx. 

In summary, we have proposed to study jet energy loss 
in high-energy heavy-ion collisions by measuring the in- 
clusive jet fragmentation function which can be extracted 
from the differential pt spectrum of charged particles in 
the opposite direction of a tagged direct photon. The 
background to the jet fragmentation is estimated to be 



small for moderately large pt. We have also demon- 
strated that modification of the jet fragmentation func- 
tion due to jet quenching can be used to obtain the en- 
ergy loss and the mean-free-path of jet interaction inside 
the dense matter produced in high-energy heavy-ion col- 
lisions. 

We have not specified the energy dependence of the 
energy loss in our calculation. In addition, the energy 
loss, dE/dx, might also depend on the distance that 
jet partons have traveled as indicated by a recent study 
[Q. These dependences can be studied experimentally 
by varying the energy of the tagged photons in the colli- 
sions of different nuclei. These are the subjects of further 
investigations 0. 
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